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1. Introduction 

Electromagnetic polarizabilities are fundamental properties of hadrons which figure promi- 
nently, e.g., in soft Compton scattering. In the low energy limit, when the wavelength of the 
scattered photons is large compared to the hadron size, the Compton scattering amplitude can be 
cast in terms of a gradient expansion in the electromagnetic field [p. The leading terms in this 
expansion are determined by the static electric and magnetic polarizabilities, which are defined 
via the hadron mass shift in constant external electric and magnetic fields, respectively. In this 
work, specifically the static electric polarizability of the neutron is investigated, building on earlier 
exploratory work Q|. By pushing the calculation down to a pion mass of 357 MeV, contact is es- 
tablished with chiral effective theory; the approach is fully dynamical, including all disconnected 
quark loop contributions. Related efforts have recently been reported by several groups, cf. [||-|5p. 

2. Expansion of the neutron two-point function 

In parallel to the Compton scattering amplitude alluded to above, the response of the neutron 
mass to an external electric field can be expanded both with respect to the magnitude of the electric 
field as well as its space-time variation. Here, only the leading term in this combined expansion is 
of interest, 

m(E)-m(0) = --E 2 + ... , (2.1) 

discarding all space-time derivatives of the electric field and terms higher than quadratic in the 
same 1 . The electric polarizability a can then be extracted by evaluating the neutron two-point 
function, which determines the neutron mass, in the presence of a constant electric field; the two- 
point function can furthermore be expanded in that field from the outset, keeping only quadratic 
terms. Consider decomposing the action into the zero external field term and the perturbation due 
to the external field, S = So + S ext . Then the two-point function expands as 

(Np(y)Na(x)) = J[DU][Dyf][D\fr] exp(-S ) {l-S ext + S 2 xt /2 + ...)N p (y)N a (x) 

J [DU] [Dy] [Df] exp(-So) (l - S ext + S 2 ext /2 + ...) (2.2) 

As usual, the terms in the denominator cancel statistically disconnected parts of diagrammatically 
disconnected contributions. The coupling to the external electric field still remains to be specified. 
In the following, the electric field will be taken to point into the spatial 3-direction, and will be 
introduced via the 3-component of the gauge field, A3 = E(t — to). This choice will be discussed 
in depth further below; it decisively influences the physics of the external field method calculation 
pursued here. For now, it is sufficient to remark that only A3 is non-zero, that it is constant in the 
3-direction, and that it is linear in E. Under these circumstances, the gauge links are modified by 
the external field as 

'From the outset, the absence of a neutron electric dipole moment is assumed in the present treatment; this would 
of course no longer be accurate in the presence of a non-trivial vacuum ©-angle. 



2 



Progress toward the chiral regime in the neutron electric polarizability 



Michael Engelhardt 



O 00 o 

II It I I I I I I I I I I 



Figure 1: Diagrammatic representation of contributions to the neutron electric polarizability. Crosses denote 
interaction vertices linear in the external field, circled crosses interaction vertices quadratic in the same. 

(where a denotes the lattice spacing and q the quark electric charge matrix) and, inserting this into 
fermion action discretizations of the Wilson type, the coupling to the external field takes the form 2 

Sext =Zv\Y< V^) t(^ A 3 - a 2 q 2 A 2 3 /2) ■ U 3 (x) ■ (-1 + fc) • y(x + e 3 )) (2.4) 

+ (iaqA 3 + a 2 q 2 A]/2) ■ U% (x - e 3 ) ■ (1 + ft) • y (* - e 3 )) 
Thus, two interaction vertices are generated, one linear in the external field, and one quadratic. Car- 



rying out the fermionic integrations in ([2.2|), i.e., applying Wick's theorem, leads to these vertices 
being inserted into quark propagators, yielding the diagrammatic representation of the relevant 
contributions to the neutron two-point function depicted in Fig. [[]. 

3. Role of the external gauge field 

As already mentioned in the previous section, the external electric field will be introduced via 
a non-vanishing 3-component of the gauge field, 

A 3 =E(t-t )=A + Et . (3.1) 

To arrive at a cogent interpretation of the neutron two-point function measured in the presence 
of such a field, a detailed consideration of the resultant physical effects is called for. First and 
foremost, note the freedom in the choice of the origin in the time direction, to, which is equivalent 
to a freedom of adding a constant, A, to the gauge field. If the spatial directions were infinite, 
this choice would be inconsequential, since A could be freely shifted by a gauge transformation. 
However, on a finite space of linear extent L with periodic boundary conditions, only a discrete 
gauge invariance A — > A + ln/qoL remains, where qo = e/3 represents the elementary electric 
charge present in the system. Apart from this residual invariance, physics, and, in particular, the 
mass shift of the neutron depends on the chosen value of A. The question therefore arises how the 
electric polarizability proper can be disentangled from the A-dependence of the neutron mass. 



2 The renormalization factor Zy is due to a technical detail: In the numerical calculation, quark lines in Fig. |l| will 
be populated by domain wall fields; however, the coupling to the external field is not realized via the five-dimensional 
conserved current, but via the current obtained by projecting the quark modes onto the four-dimensional domain walls. 



The fields iff, x// in ( 1.4 ) are the projected ones. This amounts to a renormalization of the vertices, which is compensated 
by Zy. The renormalization factor is determined to be Zy = 1.14 for light quarks and Zy = 1.10 for strange quarks by 
measuring the number of valence quarks in baryons. 
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Physically, qA enters the quark Hamiltonian 3 as a quark Bloch momentum in the periodic 
lattice. Thus, for general A, there is a current of up quarks flowing around the torus with each up 
quark, on the average, carrying momentum 2eA/3, and a current of twice as many down quarks 
with, on the average, half as much momentum each in the opposite direction 4 . Only the complete 
(periodic) "neutron", if it can be viewed as such, maintains vanishing momentum. Of course, 
the notion of a neutron becomes tenuous when the Bloch momenta become large; in this limit, 
the system essentially consists of two strong, oppositely oriented quark electrical currents with 
some residual correlations due to the chromodynamic force. Certainly, any mass shift extracted in 
such a regime could hardly be construed as representing the polarizability of a reasonably isolated 
neutron. Instead, it is necessary to seek out the regime where the quark currents cease and a bona 
fide neutron, only weakly distorted by the electric field E proper, is formed. 

In the absence of an actual electric field, it is clear how to achieve this situation; one must 
set A = 0. When E ^ 0, on the other hand, the situation is more complicated, since the Bloch 



momenta change (linearly) in time: In view of (3.1), a shift in the time direction is equivalent to a 



shift in A. This is, of course, entirely natural; the electric field accelerates the up and down quark 
electrical currents. A practical and model-independent way to identify when the quark currents 
cease is to look for stationarity in physical quantities such as the mass shift. If the quark currents 
are initially oriented opposite to the acceleration provided by the electric field, then time evolution 
will lead to a point at which the currents cease and turn around. Expanded around that point 
in time, the leading time dependence of physical quantities is quadratic. Equivalently, one may 
perform measurements of the neutron two-point function in a fixed time region, but for different 
values of A, and subsequently seek out the stationary point as a function of A. In practice, this calls 
for measurements at three different A, permitting the determination of a parabola in that variable, 
as well as its stationary point. 

Before continuing the main line of discussion by turning to the extraction of the neutron mass 
shift from its two-point function, cf. section ^ a separate concern associated with the choice ( [Q| ) 
of the external gauge field should be addressed. While this gauge field is constant, and thus periodic 
in the spatial directions, it must exhibit a discontinuity at some point in the time coordinate in order 
to be periodic in that direction as well. At the discontinuity, the electric field exhibits a spike, with 
the potential of distorting the measurement of the neutron mass shift, which is, after all, intended 
to be taken in a constant electric field. As demonstrated in such an effect can indeed be 
seen if one insists on performing the mass measurement in a time region including the spike. For 
selected magnitudes of the electric field, the distortion can be avoided, because what is physically 
relevant is not the gauge field A3 itself, but e' qAi ; thus, if one arranges the discontinuity to be an 
appropriate multiple of 271, it has no physical consequences. However, ultimately, this issue is 
moot, since nothing forces one to perform measurements near the gauge field discontinuity 5 . In 



3 In more precise terms, this would be an effective constituent quark Hamiltonian; the electrical current itself depends 
on the effective mass experienced by the charge carriers as they propagate with a given Bloch momentum. 

4 In addition to these valence currents, there are, of course, also analogous sea quark and antiquark currents. 

5 Note, however, that for the alternative choice of gauge field Aq = Ext,, the analogous objection would indeed be 
more serious; during the mass measurement, one cannot prevent the neutron from propagating to the spatial location 
of the electric field spike. Indeed, if one projects onto definite neutron momentum, all spatial positions are explicitly 
weighted equally. 
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fact, in practice, it is not even necessary to extend the support of A3 over all times; in the present 
investigation, just as in a real laboratory setting, the external field is merely switched on a little 
before the neutron source is introduced, and switched off some time after it is annihilated. This 
does not alter the result of the mass measurement, as tested in |§] ; however, it significantly reduces 
the statistical fluctuations in disconnected diagrams, since quark loops far in the past and the future 
of (and thus not correlated with) the neutron are eliminated. 



4. Extracting the mass shift from the two-point function 



As should be clear from the discussion in the previous section, the neutron Hamiltonian in the 
presence of the external field adopted here is not time-independent; a shift in time corresponds di- 
rectly to a shift in the constant A in the external field ( p.l[ ). Before addressing the time dependence, 
it is useful to review how the mass shift would enter the neutron two-point function in the case of 
a time-independent Hamiltonian depending on two arbitrary small parameters E and A. In such a 
case, projecting onto zero-momentum, unpolarized neutrons, the two-point function behaves as 

l + 7o 



G(/7 = 0,/)=£Tr 



(N(y)N(x)) 



W exp(— rar) 



(4.1) 



for large times, where both the overlap between neutron source and true ground state W and the 
neutron mass m can be expanded in E and A, 

W = W Q + W W (A,E)+W {2) (A,E) + ... (4.2) 
m + m (2) (A,£) + ... (4.3) 



m 



Inserting this into (|4JJ), the term of second order in E and A in the two-point function becomes 



G {2 \p = 0,t) ^> W exp(-m ?) 



'W®(A,E) 
W 



in 



(2) 



(A,E)t 



(4.4) 



R(t) 



implying that the neutron mass shift to second order in E and A can be obtained from (minus) the 
temporal slope of the correlator ratio R(t) = G (2) (0/G (0) (f). 

Now, the Hamiltonian of interest here is 
not time-independent; however, if it is 
sufficiently weakly time-dependent, the 
system will adjust adiabatically while 
W® and vary slowly, and ( |4.4[ ) still 
applies 6 , except that W^ and are 
now time-dependent. This implies that 
mP^ cannot be simply read off as the 
slope of R(t) anymore; this correlator ra- 
tio generically will behave as displayed 
in Fig. |[ However, as discussed further 
above, to obtain the electric polarizability 
of the neutron, one in particular needs 




Figure 2: Generic behavior of the correlator ratio 
R(t)^G (2 \t)/G^(t). Regions A, B and C qualitat- 
ively correspond to the time windows displayed in Fig. 3. 



6 A treatment beyond the adiabatic approximation, in the related context of the acceleration of a proton in an external 
electric field, has been described in \m. 
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Figure 3: Correlator ratio R(t) in (Gaussian) units of E 2 a 4 ; in these units, the plots are independent of any 
actual value of E used in the numerical computation. In the three measurements, the neutron source was 
placed at the fixed location t = 0, and the background gauge field A3 =E(t — to) was varied, with, from left 
to right, fo = 6a, ?o = and fo = —10a. This generates time windows corresponding qualitatively, from left 
to right, to regions A, B and C in Fig. 2, respectively; note the variation in the vertical scale. The behavior 
indicated in Fig. 2 is thus demonstrated, with the inflection point, corresponding to the stationary point of 
the neutron Hamiltonian, being captured in the left panel. 



only the mass shift at the stationary point in A, or, equivalently, t. Expanded around that point, the 
time dependence of and is quadratic; therefore, at that point, and only there, the mass 
shift mP^ is still given by (minus) the slope of R(t), thus permitting the extraction of the neutron 
electric polarizability. 



5. Results 

Numerical measurements were carried out in a mixed action scheme, using 448 dynamical asq- 
tad quark configurations provided by the MILC collaboration [g], with lattice spacing a = 0. 124 fm 
and quark masses ami = 0.01, am s = 0.05, corresponding to a pion mass of 357MeV. Quark lines 
depicted in Fig. [j] were populated with domain wall quarks; all disconnected contributions were 
taken into account and were estimated using bulk complex Z(2) stochastic sources. Fig. ^displays 
three different measurements of the correlator ratio R(t) = (t) / G^ ' (?) performed with different 
shifts A, thus probing different time windows of the behavior displayed in Fig. |[ As a result of 
using Wuppertal smeared neutron sources and sinks optimized with respect to their overlap with 
the true neutron ground state, as constructed in [Q], excited states are already negligible for t >4a 
(where t = denotes the neutron source location). On the other hand, statistical uncertainties are 
under reasonable control up to t = 10a. In practice, the average slope of the correlator ratio R(t) 
for 5a < t < 9a was determined in each of the three time windows displayed in Fig. ^ these three 
values define a parabola as a function of A, the extremum of which yields (minus) the mass shift at 
the stationary point, —m^ 2 \A stat ,E). The latter was corrected for curvature effects averaged over 
when determining the slope in each window. Alternative analysis schemes were explored and yield 
consistent results. Finally, supplementing [— mP-\A stat ,E)]/(E 2 a 3 ) with a factor 2 from ( j2~I[ ) and 
a factor — 1 from the continuation E — > iE to Minkowski time yields the polarizability, in Gaussian 
units, 

a = (0.9±1.6)-10- 4 fm 3 . (5.1) 

Fig. |] puts this result in relation to chiral effective theory, specifically the small scale expansion 
(SSE) expectation with which it is quantitatively compatible. Although the statistical uncer- 
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tainty in (5.1 ) is sizeable due mainly to the discon- 
nected contributions, the present calculation thus 
establishes contact with the pion mass regime in 
which chiral effective theory becomes applicable; 
an earlier exploratory study [Q] at m % = 759 MeV, 
the result of which is also included in Fig. ||, did 
not yet allow for a quantitative comparison 7 . 

It is again worth emphasizing that, to arrive 



Figure 4: Electric polarizability of the neutron; 
solid circles with error bars denote lattice 
measurements from the present work (at m z 
= 357 MeV) and from [§] (at m K = 759 MeV). 
The curve represents the SSE expectation 
constructed purely from experimental data. 
Apart from the pion mass, which is varied 
in the graph, the SSE expression contains a 
collection of further parameters, among them 
the nucleon and delta masses, f n and g^, cf. 

For these parameters, lattice measure- 
ments on the same ensemble at m K = 357 MeV 
exist; substituting those values into the SSE 
expression yields the open square. This can 
be viewed as the result of an alternative, 
partially resummed effective theory. 



at the result (5.1), it was crucial to carefully take 
into account the physics associated with an exter- 



nal gauge field of the form (|3.lD. As explained in 
detail in section ||[ on a finite spatial lattice, the 
external field induces opposite Bloch momenta for 
the u- and cf-quarks, with a bona fide neutron only 
having an opportunity to form at a specific time, 
namely the stationary point at which the Bloch cur- 
rents reverse direction under the influence of the 
electric field. That is where the electric polarizabil- 
ity can be extracted. The measurements displayed 
in Fig. [3] illustrate how, in a background field cal- 
culation of this type, the polarizability emerges as 
a rather small residual effect compared to the sub- 
stantial Bloch momentum dependence. 
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